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Abstract. Periodically- modulated potentials in the form of light fields have 
previously been applied to induce reversible phase transitions in dilute colloidal systems 
with long-range interactions. Here we investigate whether similar transitions can be 
induced in very dense systems, where interparticle contacts are important. Using 
microscopy wc show that particles in such systems are indeed strongly affected by 
modulated potentials. We discuss technical aspects relevant to generating the light- 
induced potentials and to imaging simultaneously the particles. We also consider 
what happens when the particle size is c;oniparablc with the modulation wavelength. 
The effects of selected modulation wavelengths as well as pure radiation pressure are 
illustrated. 
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1. Introduction 

One of the defining features of soft matter systems, their intermediate, or mesoscopic, 
lengthscale, places them at a remarkable confluence of physical properties, where they 
are susceptible to influence by "everyday" external forces. Not only are they "soft" in 
their response to applied shear, and their constituents sufficiently buoyant to display an 
interesting interplay between gravity and thermal (Brownian) motion [1], colloids are 
also affected by the radiation pressure-ordinarily thought of as negligible in laboratory 
experiments-arising as they scatter light which impinges on them. 

The coincidence of the magnitude of these forces (Brownian, gravitational, and that 
arising from scattered EM radiation) in these systems, set fundamentally by the relative 
magnitudes of Boltzmann's constant kb, the gravitational constant G, the speed of light 
c, and Planck's constant h, is compounded by a coincidence of the particle size with the 
wavelength of visible light. The remarkable consequence is that we are able to study, in 
real space, Brownian particles under the influence of light-induced potentials. Ashkin 
is credited with first realising this possibility [21 |3l IH |5l El [7] , leading ultimately to his 
development of optical tweezers [H] . These are now a widely-used tool in physics [U] and 
biology [ini [m [E] • Studies using optical tweezers generally consider very deep potential 
minima, where objects are tightly trapped; the present work considers in addition the 
case where potentials do not necessarily irreversibly capture particles but do bias them 
towards certain locations. 

Colloids have found favour as models of atomic systems, among their advantages 
being the possibility of imaging directly large colloids with light microscopy, the 
relatively slow speed at which they diffuse, and the relative ease with which interactions 
can be tailored [I3l [H]. In the apparently simple case of hard-sphere interactions, 
colloidal systems show the predicted entropy- driven fluid-solid transition at high density 
[15] and, at higher densities still, a glass transition instead of the expected crystalline 
equilibrium state, again in analogy with atomic and molecular glass-forming systems 
[T^ [TBI ITTf UHl HH]. Similar behaviour has been seen in other systems, for example 
charged particles interacting via long-ranged electrostatic repulsion [20] . 

Fluid-solid transitions are ubiquitous in nature, and are frequently the expected 
equilibrium state. Often, solid phases occur upon cooling, whereupon distinct 
density modulations are enhanced, resulting in translational and orientation order. 
Periodically modulated potentials can have a similar effect, giving rise to a controllable 
phase transition effected at will by application and removal of an external potential. 
Chowdhury was the first to demonstrate this, in two dimensions, in so-called laser- 
induced freezing (LIE) experiments [211 122] . They were followed by further experimental 
[231 [21 123 126] and theoretical [23 [2H1 123 130] confirmation. In LIE, a one- 
dimensionally-modulated potential gives rise to a two-dimensionally-modulated crystal. 
The registration of adjacent rows is mediated by Brownian motion perpendicular to the 
direction of the potential, which the particles continue to undergo, although they sit 
on average at the minima of the potential. Eor deeper potentials, the lateral motion 
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of the particles becomes sufficiently restricted that the system is only modulated in 
one dimension, but remains liquid-like along the direction of the fringes; laser-induced 
melting (LIM) occurs [3ll |32l |33l [HI |35l [Ml [STJ |3Hl [39] . There are many open questions 
relating to light-induced phase transitions, ranging from the kinetics of LIF and LIM, 
to as yet untested predictions for systems with more complex inter-particle potentials 
[lO] . time- varying potentials [H], and binary mixtures |12]. 

LIF and LIM both represent transitions between equilibrium states, brought about 
by the application of external, laser-induced potentials. In all of these cases, removal 
of the potential leads to a re-establishment of the initial state. In contrast to reversible 
light-induced transitions, there is the possibility that externally-applied modulated 
potentials may induce expected but unobserved transitions to ordered states. A prime 
example is the hard-sphere glass transition described above, though other glass and gel 
transitions may be relevant. Dynamical arrest due to co-operative effects such as caging 
is frequently offered as a possible explanation for the glass transition, e.g. [HI [ISj [S]. 
It is known that shear can lead to ordering in glasses and gels |l5l [IHl [ITJ [IHl 1^ . 
though the underlying mechanism is still far from understood [501 El] • (Note that shear 
can also lead to melting of crystals [S21 If the emerging picture of co-operative 

arrest is accurate, and given the above evidence of external field-induced ordering, 
it seems plausible that periodic light fields might also induce ordering in very dense 
systems. Biroli et al. have developed an inhomogeneous mode-coupling theory (MCT) 
which predicts the response of a supercooled liquid's dynamical structure factor when 
exposed to a static inhomogeneous potential |54j- Agreement of these calculations with 
experimental results would represent a powerful test for MCT. 

This paper describes a ffist investigation into whether periodically modulated light 
fields can indeed induce an effect in very dense colloidal systems. In a few initial steps 
towards the grand goals mentioned above, we present some technical aspects relating to 
the experimental realisation of an apparatus which permits exposure of a sample to the 
light fields and its simultaneous observation. Our first results with this apparatus show 
that even the most dense samples do show a clear rearrangement under the influence of 
an applied light field. 

2. Light as an external potential for colloidal particles 

Light forces acting on small particles have been described in considerable detail, both 
from a ray-optics and a Rayleigh perspective. The necessary basic physics is well 
established [551 EHl ISZl IS^- We outline the two features most important to the 
present work, the so-called scattering force -Fscat, which acts in the direction of beam 
propagation, and the gradient force Fgj-ad, which depends on the shape of the light 
intensity distribution. 
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2.1. Scattering force Fscat 

Historically, the scattering force preceded full optical trapping in the form of Ashkin's 
levitation experiments [21 El H] . It arises as a result of momentum transfer to a particle 
from incident photons as they are scattered from it. 

A photon of light with frequency z/, wavelength A and phase velocity vp carries 
energy -Ephot = hi/ = hvp/X and momentum p = h/X = Ep^ot/vp. In absorbing a 
photon, an object therefore experiences a force Fgcat = dp/dt = d {Ep]^ot/vp) /dt = 
P/vp = nP/c, where P is the power of the photon source. The last equality follows 
since the phase velocity of light in a material of refractive index n, vp = c/n, where 
c is the speed of light in vacuo [HI [9j|} In our case, photons are scattered rather than 
absorbed and therefore transfer only a portion of their momentum to the particle; in 
the case of dielectric spheres, a prefactor g ~ 0.1 is customarily assumed [21 E]- 

In our experiments, we have calculated -Fscat to be on the order of a tenth of the 
particles' own weight for each mW of applied laser power (as measured at the laser, 
as always in this article); in a typical experiment this corresponds to around a tenfold 
increase in the particles' effective buoyant mass. 

2.2. Gradient force Fgj-s,d 

The gradient force -Fgrad allows a single beam, suitably focused, to act as an optical 
trap. Its origin is intuitively clear in the ray optics formulation [HI EH], although since we 
ultimately consider the sum of forces over all infinitesimal volume elements of finite-sized 
spheres, we discuss the origin of the gradient force as it applies to Rayleigh particles. 
In this regime, the electric field over a particle is approximately uniform, and induces a 
dipole moment p = aE, where a is the polarisability of the (dielectric) particle given 
by the Clausius-Mosotti relation [60] : 



with n = ric/ns the ratio of the refractive indices of the colloidal particle, n^, to the 
surrounding medium (solvent), Ug, and a is the particle's radius. Having acquired a 
dipole moment, the particle experiences a Lorentz force fL = (p ■ V)E + {l/c)dp/dt x B. 
This can be re- written 1551 : 



The second term is the scattering force Fs^at discussed above, and is directed along 
the direction of propagation. The first term, Fg^ad = (a/2)VF^, is the gradient force 
and expresses that in a non-uniform electric field, a particle will move towards regions 
of higher electric field if nc>ns, and vice versa. This indicates that the modulated 
potentials we seek can be realised by spatially-modulated electric fields. 

I We have assumed the Minkowski form, which differs from the Abraham form by a factor n^. Ahhough 
this is an important unresolved theoretical issue this subtlety does not concern us unduly, since 
optical forces can be calibrated experimentally (and n^^l). 
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2. 3. Modulated potentials from modulated light fields 



One conceptually simple means of generating modulated light fields is using a crossed 
beam experiment, as employed in previous studies similar to ours [22|l33]. and in thermal 
diffusion forced Rayleigh scattering studies [611 162]. The experimental arrangement is 
also similar to the so-called "dual-scatter" or "dual-beam" configurations used in Laser 
Doppler Anemometry [631 |6l]- For two crossed coherent TEMqo mode (Gaussian profile) 
laser beams, the resulting intensity profile is [21]: 

I{x) = 2/o {1 + cos [2kxsm{e/2)]} e-2-'-«'W2)/fl^ (1) 

where 6 is the beam crossing angle, Jq the intensity of each beam, R the laser beam 
radius (the e"^ point), and k = 2tt/X the incident beam wavevector. The term in 
braces, {1 + cos [2A;x sin(6'/2)]} = {1 + cos [qx]}, is the sinusoidally-varying interference 
pattern of fringe spacing d = X/ (2sin(^/2)) and fringe wavevector q = 27i/d. The last 
term represents the overlying Gaussian envelope due to the finite beam size. Figure [Tj 
compares a calculated fringe pattern (a) with an experimentally observed pattern (b). 
From the gradient, the magnitude of the force can be calculated i.e. Fgrad = («/2)V£^^ 
(c). These images illustrate that not only are particles drawn to the fringes, but they 
are also confined by the Gaussian envelope. Figure [2] illustrates this in the form of a 
vector field plot, which was obtained by numerical differentiation (Appendix 1). 




Figure 1. Fringe pattern calculated for typical experimental parameters which result 
in a fringe spacing d = 9.29 fim (a), an observed fringe pattern (b), and the magnitude 
of the gradient force calculated based on the fringe pattern shown on the left (c). 
Parameters are (Section 3.1): wavelength A = 532nm, crossing angle 9 = 3.28°, beam 
radius at the sample R = 1.22\f/D = 23.2/im [21 for a laser beam of diameter 
D = 2.80 mm focused by a lens of focal length / = 200 mm. 



From the first term of the modulated light field. Equation [T| we obtain the 
oscillatory part of the potential: V{x) = Vq { 1 + cos [ga;] } , where Vq absorbs the 
particles' polarisability and the laser beam intensity. Loudiyi et al. additionally 
normalise this quantity by the thermal energy [21]. 
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Figure 2. A vector plot showing the position-dependent force experienced by a 
particle, found by numerical differentiation (Appendix 1) of the calculated pattern 



in Figure [T] (left). 



2.4- Effect of finite particle size 

Implicit in all of the previous section was that the potential acts only at the centre of 
the particles, or alternatively that the particle radius a <^ d. This is not satisfied in the 
present work. For finite particle size, integration over infinitesimal volume elements of 
the particle results in a modified potential |24j : 



where x is the position of the particle centre and ji the first order spherical Bessel 
function. 

The finite size of the particles leads to the additional factor 3ji(ga)/(ga) (Figure 
[3]). This factor resembles the form factor of a sphere as it occurs in scattering 
experiments [65]. In a scattering experiment, the sample is illuminated by a beam and 
the intensity scattered under a scattering angle 9 or, equivalently, a scattering vector q 
is determined. This implies an identical geometry with the incident and scattered beam 
here represented by the two crossing beams. 

The additional factor demonstrates several interesting features [66]. First, for large 
fringe spacing [qa ~ 0), the particles behave like point particles. Conversely, for very 
small fringe spacing (qa^oo), the effective potential is averaged to zero. Interestingly, 
there are fringe spacings where the sign of the potential is reversed, indicating that 
spherical particles can be either drawn into the fringes or repelled from them, depending 




(2) 
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Figure 3. The factor 3ji{qa)/{qa) arising from the finite particle radius a (Equation 



on qa. As long as the fringe spacing is greater than d = 0.699 x 2a (corresponding to 
qa = 4.493, the first root of ji{qa)), the factor 3ji(ga)/(ga) > 0. This means that 
for fringe spacings larger than the particle diameter but still finite (and indeed slightly 
smaller too), the behaviour is qualitatively similar to that for point particles, albeit with 
reduced potency. 

3. Experiment 

3.1. Apparatus 

Figure |4] shows the experimental arrangement used to create the modulated potential 
and, at the same time, observe the response of the sample. The modulated potential is 
created by splitting a linearly-polarised laser beam (Coherent Verdi V5 with P = 5 W, 
A = 532 nm) and subsequently crossing the two beams. The beam is split using a 50:50 
beamsplitter (BS) with a preceding half- wave plate (A/2) to adjust the polarisation for 
optimum performance of the beamsplitter. The two beams are brought parallel to one 
another by means of two mirrors (Ml, M2) and a moveable pair of mirrors. Translation 
of the moveable mirrors adjusts the beam separation s and, after the focusing lens (LI), 
the crossing angle 6 of the two beams. A half- wave plate (A/2) in one of the beams allows 
rotation of the polarisation of one beam with respect to the other, thereby controlling 
the amplitude of the interference fringes whilst maintaining a constant mean intensity 
and thus radiation pressure. 

The introduction of the sample changes the crossing angle. For a typical sample 
cell (depth 170 /im) and fringe spacing (7 /im), the angle of incidence is reduced from 
ei/2 = sin"^ {X/2d) = 2.2° to er/2 = sin"^ (sin(2.2°)/1.33) = 1.6°, corresponding to 
change in the focal position of 170 fim x (tan(^j/2)/tan(^r/2) — 1) ~ 64/im. This is 
corrected by a linear translation of the lens LI (Figure |4]). Note that despite the change 
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in the crossing angle, the fringe spacing remains unchanged (since A also changes upon 
entering the new medium.) 



LAMP 




OPTICAL FILTERS 

Figure 4. Schematic diagram of the experimental arrangement used to generate a 
modulated light field and simultaneously image the sample. Two coherent beams are 
created by beamsplitter BS, and subsequently brought parallel by two mirrors (Ml, 
M2) and a pair of moveable mirrors. The position of the moveable mirrors determines 
the beam separation s, and therefore the beam crossing angle 9 following focusing lens 
LI. For observation, white illuminating light is combined from above using dichroic 
mirror Dl and the sample imaged by an objective. The majority of the intense laser 
light is deflected to beam dump BD using a second dichroic mirror D2. A series of 
three optical filters, some rotatable, in the objective infinity space allows for variable 
attenuation of the remaining laser light and thus for an adjustment of the brightness 
of the fringes to be viewed. Images are recorded, after suitable magnification, using a 
digital camera. 

Concurrently, the sample is imaged with a home-built inverted bright-field 
microscope. Kohler illumination is provided from above the sample. The extremely long 
working distance condensing lens provides sufficient space for a "notch" dichroic mirror 
(Dl), reflective in a narrow range around A = 532 nm but otherwise transmitting in the 
visible. This dichroic mirror combines the imaging and modulated light at the sample. 
After the sample, a standard high numerical aperture microscope objective (Nikon xlOO 
PA VC, NA=1.4) forms an image at infinity, before a tube lens (TL) and subsequent 
telescope (L2, CO) adjust the magnification as appropriate for the camera. Additional 
optics can be introduced straightforwardly into the (relatively long) so-called infinity 
space behind the objective. We use this to separate the intense laser light, damaging 
to the camera, from the white imaging light. The bulk (~ 98%) is deflected to a beam 
dump (BD) using a second dichroic mirror (D2). As well as eliminating safely the 
majority of the laser light, this light can be re-used, for example by retro-refiection to 
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achieve a counter-propagating arrangement whereby Fgcat can be reduced independently 
of the laser intensity. Even after the second dichroic mirror, the intensity of the laser 
light is far too high for the camera. Three additional optical filters are used to adjust 
the level of the modulated light field while retaining most of the imaging light. This 
permits simultaneous imaging of the sample and an appropriate fraction of the light 
field. 




Figure 5. Variation of the filter optical density at 532 nm, O-D532, with the angle 
of incidence, 7. Pq and P are the measured power values before and after the filters, 
respectively. 

The intense laser light has to be reduced to about 10"^ - 10~* W at the camera, to 
obtain satisfactory images of the fringes. The fringes become essentially invisible upon 
a further reduction by a factor of about 100. Since the laser power varies depending 
on the experiment, for optimum simultaneous imaging, filters with a variable optical 
density at 532 nm, OD532, are desirable. This can be achieved for interference filters 
by changing the angle of incidence, 7 (Figure |4]): the filters have a sharp transmittance 
edge slightly above 532 nm which shifts to lower wavelength as 7 is increased (similar to 
notch filters [67]). Since the slope of the edge is finite, this provides control over OD532 
(Figure [5]). When placed in the infinity space of the microscope, these filters can be 
straightforwardly rotated to allow the fringes to be imaged or not, as desired. The weak 
dependence of O-D532 within — 5° < 7 < 5° (Figure [s]) is important for imaging, since 
this range is slightly larger than the divergence in the infinity space of the microscope 
[§} Nevertheless, for different parts of the field of view, the effective 7 and thus OD532 
is different, and hence an image of the fringes is no longer quantitatively correct. A 
correct image can, however, be obtained with 7 = 0°, or by using neutral density filters 

§ Light from the focal plane is focused at infinity, but, except for light originating from the point on 
the optical axis, is nonetheless divergent. The focal length of an objective is the microscope tube lens 
focal length (here about 200 mm) divided by its magnification (here 100). Together with the radius of 
the field of view (here about 125 ^m), this results in a divergence in the infinity space of the microscope 
of around tan~^(125 /im/2 mm) ~ 3.6°. 
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(in which case the white hght is attenuated beyond usefulness). The bright-field images 
remain good since the transmittance of each filter at wavelengths A 7^ 532 nm is T ~ 0.9. 



3.2. Analysis 

Having determined the particle coordinates [68t [69] . a range of parameters can be 
calculated, for example the particle density 0, the pair correlation function g{r), the 
mean coordination number (z), the distribution of coordination numbers p{z), and 
bond-orient at ional order parameters, e.g. Calibration of distances, necessary 

for determining the fringe spacing as well as for structural analyses, is performed using 
a high- resolution microscope test slide (Richardson Test Slide, Model 80303) [69] . 



3.3. Samples 

We have used polystyrene sulphate spheres of radius a = 2 fim (Interfacial Dynamics 
Corporation) suspended in water. The large refractive index difference between particles 
(ric = 1.59) and water (ug = 1.33) results in large optical gradient forces, but the 
concomitant multiple scattering limits their use effectively to a single layer, i.e. two 
dimensions. These spheres carry negative charges, which in the present study are 
screened by high salt concentrations. We regard them as (almost) hard spheres, which 
is supported by the observed distance of closest approach and the shape of the pair 



correlation function (Section 5.2). Though the salt concentration is high, it is still low 
enough to avoid problems with coagulation. 

Samples are prepared by pipetting a suitably-diluted homogenised stock solution 
directly into the sample cell, which fills largely by capillary action. The particles quickly 
sediment on the coverglass. For dilute samples, this results in two-dimensional samples, 
while at higher concentrations a few layers form (which can be reduced to a single layer 
by application of radiation pressure. Section |4]). The concentrations we refer to in the 
following are the volume fractions of the initial, homogenised, bulk solutions. This is a 
nominal value; in the final sample, inhomogeneities in the density may occur depending 
on the settling process. 

The sample cells consist of coverslips glued together, giving a sample volume of 
about 20 mm x 3 mm x 170 /im ||69j. The sample only comes into contact with glass 
and possibly the UV-cure glue used to seal the cells, whose effect is assumed negligible. 
Since glass becomes negatively charged in the presence of water, the particles are repelled 
from the surfaces of the cell, and become attached only very occasionally [TTj [72]. 



4. Radiation pressure results 

We first investigate the effect of radiation pressure only as a function of laser intensity 
and particle concentration (Figure |6|. The radiation pressure is applied by turning the 
second half- wave plate (A/2) until minimum contrast is achieved as judged from images 
of the interference patterns formed using neutral density filters. The concentrations are 
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chosen such that a first layer of particles next to the coverslip (slightly out of focus in 
Figure |6]) as well as an incomplete second layer (in focus) are formed. 



Figure 6. The effect of pure radiation pressure as a function of sample concentration 
(from top to bottom: initial homogenised volume fractions cf) — 0.015, 0.020, and 0.030) 
and laser intensity (left to right: P=0.00 W, 0.10 W, 0.20 W, and 0.50 W). Images 
are taken following 20 min of irradiation, except for the top right image which shows 
the sample (0 = 0.015, P=0.10 W) 30 min after the laser is turned off. Note that only 
the particles in the second layer are in focus; they show a bright spot at their centre. 
Particles in the first layer (most prominent in the image second from left, top row) are 
out of focus but still clearly identifiable. 

As the sample concentration increases (downwards in Figure |6]), the density of the 
second layer increases. The presence of a second layer does not imply that the maximum 
possible density has been achieved in the first layer. Indeed, this is observed not to be 
the case: upon increasing the radiation pressure, particles are pushed from the second 
layer into the first layer (left to right). In each of these cases the radiation pressure was 
applied for 20 min. In the least dense sample (top row), a laser intensity P =0.10 W 
is already sufficient to insert all of the particles into the first layer. (The two highest 
laser intensity results for the lowest concentration are omitted in Figure |6}) This forms 
a dense hexagonally- close packed (HCP) layer. With increasing density, a greater laser 
intensity is required to insert all of the particles within the laser beam into the first 
layer. Beyond a certain density, it is no longer possible to insert all of the particles into 
the first layer, even for very large radiation pressures. With increasing concentration, 
the area with only a single layer gets smaller, corresponding to the Gaussian profile of 
the laser beam and thus the applied radiation pressure. This can also be seen with 
increasing laser intensity. 
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The top rightmost image shows the least concentrated sample 30 min after the field 
is removed. (Similar relaxation behaviour is observed in all samples.) The second layer 
has become repopulated. This indicates that the osmotic pressure experienced within 
the highly concentrated first layer is sufficient to cause particles to "pop up" into the 
second layer. This also implies that the inability of particles to enter the first layer 
under their own weight cannot be explained by pure geometrical frustration. We have 
observed that "popping up" occurs with a characteristic time of about 10 s. If this 
upward movement into the second layer is thermally driven, i.e. is a chance Brownian 
excursion (and a return to the first layer is hindered by particle rearrangements within 
the first layer), its timescale should be given by Kramer's escape time with a ramp 
potential of depth Uq and extent 2a representing the gravitational potential [HI 173] : 

' dx'e^^^^') r' dxe-^^^^^ '''''' 



where /3 = l/k^T and is the self-diffusion constant for a free particle. Using 
appropriate parameters gives r~10^s^lOs, indicating the presence of a driving 
force, namely the osmotic pressure. 

This seems plausible in light of the established connection between the statistical 
geometry of hard spheres and their thermodynamic properties [TU [751 ES]. These 
references suggest that insertion of particles into a disordered layer by the application 
of radiation pressure should, as in the Widom insertion method, permit study of the 
thermodynamic properties of particles in the ffrst layer. 



5. Modulated potential results 



Having established a two-dimensional sample, we now introduce a modulated potential. 
In the present study, the modulation is always as great as possible at the speciffed laser 
intensity, i.e. both beams have the same direction of polarisation. In each experiment, 
we prepare a dense hexagonally-close packed (HCP) layer using radiation pressure (at 
the indicated power) before the half-wave plate is rotated to "turn on" the modulation. 
The parameters we vary are the fringe separation and the amplitude of the modulation. 



The HCP symmetry suggests a few fringe spacings d (Figure 12 A.); here we investigate 
d = y/Sa (Section 5.1) and d = 2\/3a (Section 5.2) 



5.1. Natural fringe spacing 

For a modulated potential of fringe spacing d = VS a, it is possible for all particles 
forming an HCP layer to lie at the potential minimum. We thus consider this a natural 
fringe spacing. 

Pure radiation pressure (a single beam of P =0.50 W for 30 min) leads to randomly- 
oriented crystallites (Figure [7| left). After exposure to the modulated potential 
(P =0.50 W for 30 min), the crystallites have rotated and consolidated to a near-perfect 
crystal with a clear direction aligned with the fringes (right). This is also reffected in the 
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Figure 7. Sample (initial concentration 4> — 0.020) after 30 min of pure radiation 
pressure (a) followed by a further 30 min with a modulated potential of wavelength 
d = ^/Sa (b). The laser intensity in both cases is P =0.5 W, and the superimposed 
line indicates the approximate fringe direction. 

distribution of the nearest-neighbour bond direction, which shows three strong peaks 
separated by 60° (Figure |8|. 

120 

100 

§ 80 
c 
B 

3 60 
u 

o 

40 
20 



"0 30 60 90 120 150 180 

Angle / ° 

Figure 8. Distribution of the nearest-neighbour bond direction for a sample (initial 
concentration 4> = 0.020) after 30 min of pure radiation pressure (solid black line) 
followed by a further 30 min with a modulated potential of wavelength d = y/3a 
(dashed red line). The laser intensity is P —0.5 W. 

The crystaUites thus seem to be able to rearrange despite the high density. It is 
interesting to investigate exactly how this process occurs. One observation is that as part 
of a crystallite rotates, the total energy in the light field does not decrease monotonically 
until the particles are aligned with the field. At some angles ip between the light field 
and the crystal orientation, relatively many particles are near to the potential minima. 
When the particles are aligned with the field {il!=0°), all of the particle centres occupy 
a minimum. For angles '?/'~18°, 31°, and 42°, there are only around 40%, 60%, and 
40% of the particles in the minimum respectively, whereas for in-between angles there 
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are far fewer. Supposing a large crystallite were to rotate towards the global minimum, 
therefore, it may do so at varying speed, perhaps even pausing at these intermediate 
metastable orientations depending on the amplitude of the field. 

5. 2. Twice natural fringe spacing 

We now consider a fringe spacing d = 2^/3 a corresponding to twice the spacing between 
two rows of an HCP layer. For sufficiently dilute samples, the particles align along the 
fringes (Figure |9]). 




Figure 9. A dilute sample exposed to a modulated light field (P=0.50W) with spacing 
d = 2^/3a. 



In dense samples, more complex structures develop (Figure 10). The initially 
disordered sample (top left) develops randomly-oriented crystallites following the 
application of radiation pressure (1 h of P =0.40 W, top right), as described previously 
(Section |4|. Relatively soon after a modulation of wavelength d = 2y/3a is introduced 
(100 s, still with P =0.40 W, bottom left), the sample is altered, with the emergence of 
voids which run broadly in the direction of the fringes. After substantially more time 
(about 5 h, bottom right), the field has caused significant structural rearrangement. 
In time-lapse movies of images, groups of clusters can be seen moving co-operatively, 
leading to arrangements along the potential minima. In particular, the motif highlighted 
in Figure 10 (bottom right) occurs frequently, with an orientation relative to the fringes 
as indicated in Figure [T2| d. This rotation is understandable on energetic grounds, which 
we discuss further below. Other samples show similar behaviour. 

The structural evolution of the sample has been investigated more quantitatively 
by following the rearrangements induced by a modulated potential (Figure 11). We 
determined the positions of particles which were located in a rectangular region within 
the single-layer region and thus under the infiuence of the modulated potential. Over the 
course of the whole period, the number of particles N{t) within the observation region 
and thus the particle density steadily decreases (a). In addition, a particle's average 
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Figure 10. Micrographs of a sample (initial concentration (j) = 0.020) before 
irradiation (a), following 1 h of radiation pressure at a laser intensity of P =0.40 W 
(b), and 100 s (c) and about 5 h (d) after the introduction of fringes with spacing 
d = 2\/3a. 



number of neighbours {z{t)) drops from around 4.4 to 3.9 after 3.5 h (b). This is also 
reflected in the distribution of the number of neighbours p{z,t) (c), which indicates 
an increasing probability of weakly connected particles, consistent with the appearance 
of voids along the fringes. Although the number of neighbours decreases, the bond- 
orientational order parameter ipQ (Section |3.2 ) does not change significantly over the 
course of the experiment (d). This indicates that those particles which remain bonded 
do so in a morphologically similar way. This is supported by the fact that the pair 
correlation function g{r) is also essentially unaffected throughout the experiment (e). 

How can we understand these observations? In dilute samples, all of the particles 
can be arranged in the potential minima. For the dense samples, half of the particles 



can still lie along the minima (Figure [l2p.) (Section 5.1), but the remaining particles 
are forced to lie between the fringes and thus at the maximum of the potential. While 
the intensity gradient and hence the force is zero at the maximum, this arrangement 
is metastable, with very small fluctuations inevitably resulting in large gradient forces. 
These forces attempt to insert particles into the minima, i.e. the fringes (dark blue 
particles in Figure I12]d), and in so doing push other particles along the fringes (blue 
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arrows). This is achieved without penalty, provided the density at the end of the 
fringe is suitably low. When the density of the sample is large, there is a significant 
osmotic penalty associated with pushing particles along the fringes and into the bulk. 
A balance must be struck between the optical gradient force and the osmotic force, 
which are opposed in their preference for density modulations. This explains why the 
expected modulations in density are observed at low concentrations (Figure |9|, but not 
at very high concentration. At high densities, the system aims to accommodate as many 
particles as possible within the fringe, without significant extension along the fringe. We 
have observed structures which achieve this; one example is that highlighted in Figure 



10 (bottom right) and explained in Figure 12 i. This rhombic "motif" represents a part 



of the crystal which, after rotation through 30°, reaches an energetically advantageous 
state (which depends on the precise details of the potential, see below) without a large 
extension along the fringe direction. The rearrangement of small crystalline parts leaves 
bond orientations unchanged, consistent with the observed essentially constant ip6: cis 
well as leaving inter-particle distances largely unchanged. This latter observation is 
consistent with our finding that g{r) does not change substantially. What modest 
extension along the fringes there is expels some particles, in agreement with the decrease 
in the particle number N{t) and in turn the mean number of neighbours {z{t)). 

These observations might have interesting consequences. First, if what we observe 
are equilibrium structures, it is remarkable that they form via small crystalline parts 
which are broken away and simply reoriented with respect to the applied potential. It 
is, however, also conceivable that, due to the geometrical frustration in a dense system, 
these co-operative motions are the only means by which the system can rearrange. In this 
case, the observed structures would correspond to non-equilibrium states liable to further 
evolution; indeed the evolution of particle number N{t) and mean coordination number 
{z(t)) suggests that the samples are still evolving (Figure lli,b). Whether equilibrium 
or not, it is clear that the modulated potential has a profound effect even in these dense 
samples. Our experiments also suggest that at intermediate (in the present context, 
though these are still relatively very dense samples) concentrations, novel structures 
might form due to the competition between the imposed potential which favours density 
modulations, and the osmotic pressure of the system which opposes them. 

Which structure is energetically or kinetically preferable depends on the shape 
of the potential. For example, for a square-well potential, the particles can, to some 
extent, move laterally within the fringe without penalty. Depending on the potential 
width and separation, a wealth of structures has been predicted [77j. Although in 
that study the colloids remain near to one another due to mutual attraction rather 
than osmotic pressure (as in our case), the effect is seemingly similar. For a potential 
with monotonically increasing curvature, e.g. a quadratic potential, it is advantageous to 
displace particles from the minimum as little as possible; "zig-zag" lines are expected. In 
the present case, however, the curvature of the potential is non-monotonic and it seems 
reasonable that some particles maintain their position while others are significantly 
displaced from the potential minimum. Together with the influence of the osmotic 
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pressure due to the bulk sample, this energy-minimisation argument justifies the 
existence of the observed motifs. 

6. Conclusion 

We have described an apparatus used to expose a sample to sinusoidally-varying light 
fields and simultaneously image the sample. To demonstrate its capabilities, we have 
investigated the response of colloidal particles to the modulated potentials which arise 
from the light field. We have shown that these potentials infiuence even samples dense 
enough that the dynamics of their constituent particles are severely restricted. Currently 
we are further improving the apparatus by including a counter-propagating beam which 
will allow us to control the modulated potential and radiation pressure independently. 
This will be achieved by replacing the beam dump by a retro-refiector. 

Densely-packed effectively two-dimensional samples have been generated using 
radiation pressures of different intensity. The behaviour of these samples upon exposure 
to modulated potentials has been investigated for two different modulation wavelengths. 
This has revealed co-operative structural rearrangements and final structures which seem 
to result from a competition between the optical gradient force and the osmotic pressure 
of the bulk sample. While theoretical predictions for a sinusoidal potential are lacking, 
similar theoretical calculations suggest structures comparable to those we have observed. 

With this apparatus, we can now investigate different situations: first, for disorder- 
to-order transitions, specific predictions exist for binary hard disc mixtures under similar 
conditions to those described here Second, disorder-to-disorder transitions are 

expected for systems with attractive interactions exposed to modulated potentials |10] . 
Both these transitions represent reversible transitions starting from equilibrium states, in 
which the initial states are recovered on removal of the modulated potential. In contrast, 
in a third situation, high-density non-equilibrium systems, in particular repulsive and 
attractive glasses, might undergo irreversible transitions from their non-equilibrium state 
to an ordered equilibrium state upon exposure to a modulated potential. In this case, 
structural rearrangements lead to stable configurations that persist even after removal 
of the external potential. In addition to revealing new physics, this might also have 
implications for material sciences. 
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Appendix 1: Numerical differentiation of patterns 

The calculated intensity profile was differentiated numerically using the Sobel method 
[78| §7.1.3] to obtain an approximation to the force field experienced by the particles. 
If / is the image, then the gradient of the image 



Vf = 







df/dx 


Gy 




df/dy 



with magnitude |Vf| = {Gl + GyY^^ and direction ip{x,y) = tan ^{Gy/Gx) is formed 
by convolution of the image with the following kernels: 
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Figure 11. Effect of a modulated potential on the evolution of different parameters. 
Shown are the time dependence of (a) the number of particles N{t) within the 
observation region, (b) the average number of neighbours {z{t))^ (c) the distribution of 
the number of neighbours p{z^t) with time (direction of increasing time indicated by 
arrow), (d) the bond-orientational order parameter ijjQ{t), and (e) the pair correlation 
function _g(r, t). 
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Figure 12. (a) Hexagonally-close packed (HCP) layer of particles with radius a and 
inter-layer spacing i/Sa. (b) Particles located along the fringes (minima) are stable 
while those at the maxima (indicated in dark blue) are metastable and, as a result 
of fluctuations, experience a force toward the fringes. They can join a minimum if 
particles which are already present in the minimum can advance along it (blue arrows), 
(c) When this is hindered, the particles can locally rearrange (e.g., rotate by 30°) to 
adopt more favourable structures. 



